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Erfindungsmeldung 



Schaltungsanordnung fur einen Zeitbereichsentzerrer bei der dis 
kreten Multitonmodulation 



1. Welches techmsche Problem soil durch die Erfindung geldst werden? 

Die diskrete Multitonmodulation (DMT) kann zur Obertragung von Daten Qber linear verzer- 
rende Kanaie herangezogen werden. Dabel kdnnen sich bezGgllch Flexibility bei der Anpas- 
sung der Datenrate oder des Sendespektrums an den Obertragungskanal oder die Stfirum- 
gebung Vorteile gegenuber Eintragerverfahren ergeben. Das Verfahren wlrd fQr den digitalen 
Rundfunk (DAB = Digital Audio Broadcast) unter der Bezeichnung OFDM (Orthogonal Fre- 
quency Devision Multiplex) eingesetzt. FQr den digitalen asymmetrischen Teilnehmeran- 
schlufi Dber normale Telefonlertungen (ADSL = Asymmetric Digital Subscriber Une) wurde 
dieses Verfahren ebenfalls standardisiert. 

Bei dem Modulationsverfahren setzt sich das Sendesignal aus vielen sinusfdrmigen Signalen 
zusammen, wobel jedes einzelne Sinussignal sowohl in der Amplitude als auch in der Phase 
moduliert wird. Man erhSIt somit eine Vielzahl von quadraturamplitudenmodulierten Signa- 
len. 

Zur Reallsierung venwendet man im Sender die inverse FFT (Fast Fourier Transformation) 
und im EmpfSnger die normale FFT, fQr die es effiziente Signalverarbeitungsalgorithmen 
gibt. 

Das Bild 1 zeigt das prinzipielle Blockschaltbild des DMT-Obertragungssystems. 




Rauschen 



Bild 1: Blockschaltbild des DMT-Obertragungssystems 
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Die zu Qbertragenden Daten werden zunachst in Blocke zusammengefaftt Je nach Stufig- 
keit werden dann Jewells eine gewisse Anzahl von zu Qbertragenden Bits einer komplexen 
Zahl zugeordnet. N solcher komplexer Datenvektoren werden dann anschlleOend m'rt Hilfe 
der inversen diskreten Fouriertransformatlon (IDFT) in den Zeitberelch transformlert Die 
nach der Transformation erhaltenen Werte stellen unmlttelbar die Abtastwerte des Sendesi- 
gnals dar. Die komplexen Datenvektoren konnen somit als Amplituden von innerhalb eines 
Blocks auszusendenden Kosinus - (Realteil) und Sinusschwingungen (Imagin&rteiQ aufge- 
faBt werden. Die Frequenzen sind dabei aquidistant vertellt und betragen 

/ /t= /.l 1-1.2 N/2 

wobei T die Zeitdauer eines Blockes und N die Anzahl der Abtastwerte eines DFT-Blockes 
darstellt Wahlt man for N eine Zweierpotenz, dann kann fur die Transformation in den Ze'rt- 
bereich die inverse FFT (Fast Fourier Transform) angewendet werden, wobei der Realisie- 
rungsaufwand stark reduziert wird. *. _x «, ,, «/,^ M \ 

Vor Weitergabe an den D/A-Wandler werden die letzten M Abtastwerte eines Blockes (M<N) 
an den Blockanfang nochmals angehangt (Cyclic Prefix). Hiermit kann im Empfanger em 
periodisches Signal vorgetauscht werden, wenn der durch den Obertragungskanal verur- 

sachte Einschwingvorgang nach M Abtastwerten, entsprechend einer Zeit von T~, ab- 

geklungen ist Der Entzerrungsaufwand ISlit sich somlt Im Empfanger stark reduzleren. da 
nach der Demodulation im Empfanger nur mit dem inversen Kanalfrequenzgang muWpMziert 
werden mufi. urn die linearen Verzerrungen des Obertragungskanals zu beseltigen ( je Fre-. 
quenzlinie eine komplexe bzw. vier reelle Multiplikationen). 

Bei der DatenQbertragung Qber Kupferleitungen treten relativ lange Einschwingvorgange auf. 
Andererseits ist anzustreben, den Zyklic Prefix mSglichst gering im Verhalthis zur BlocklSnge 
zu halten ( M«N ), urn eine moglichst geringe Reduzierung der Obertragungskapazitat zu 
erhalten Zum Beisplel wurde beim ADSL-Standard for die DatenQbertragung von der Ver- 
mittlung'zum Teilnehmer eine Blocklange von N=512 und den Zyklic Prefix eln Wert von 
M=32 festgelegt Da der Einschwingvorgang nach 32 Abtastwerten noch nlcht abgeklungen 
ist treten im Empfanger zusatzliche Fehler auf. die durch den Frequenzberelchsentzerrer 
nicht kompensiert werden konnen? Das Bild 2 zeigt das mit Hilfe einer Simulation gewonne- 
ne S/N-Verhaltnis in dB nach der Demodulation und einer idealen Frequnzbereichsentzer- 
rung unmittelbar vor dem Entscheider fur alle genutzten Frequenzkanale fQr ein ADSL- 
System Von den N=256 Frequenzlinien werden die Linien 7 bis 250 fur die DatenQbertra- 
qung genutzt. Der Linienabstand betragt 4.3125 kHz. so dafi sich fQr das Obertragungfre- 
quenzband der Bereich von etwa 30 kHz bis 1.1 MHz ergibt. Als Obertragungsleitung wurde 
eine Zweidrahtleitung mit 3 km Lange und einem Durchmesser von 0.4 mm gewShlt 
(Parameter gemaB ETSI DTR/TM-06001). 

Ohne BerQcksichtigung von aulieren und inneren Starungen (Kanalrauschen. Nebenspre- 
chen, Realisierungsungenauigkeiten des A/D und D/A-Wandlers sowle des Sendeverstar- 
kers, Quantisierungsrauschen ) ergeben sich fur das S/N-Verhaltnis am Entscheider- 
Eingang Werte im Bereich von etwa -21 bis -38 dB. 

Der durch Nichtabklingen des Einschwinganteils verursachte Fehler nach der Demodulation 
kann mit Hilfe besonderer SignalverarbeitungsmaBnahmen reduziert bzw. bese'rtigt werden. 
Die vorliegende Erfindungsmeldung beschreibt ein Verfahren und eine Schaltungsanordnung 
zur Beseitigung derartiger Storungen. Der Realisierungsaufwand gemessen in Anzahl der 
Multiplikationen je Sekunde fur die Korrektur sowie fur die adaptive Einstelluhg ist bei die- 
sem Verfahren erheblich geringer als bei den z.Z. verwendeten Verfahren . 
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Bild.2: Signalstorabstand am Eingang des Entscheiders for atle genutzten Frequenzkan§le, 
ohne Zeitbereichsentzerrer, mit idealem Frequenzberelchsentzerrer, 
Leitungsiange 3 km (0.4 mm Durchmesser) 



2. Wie wurde dieses Problem bisher gelost ? 

Im Bild 3 ist das Blockschaltbild eines DMT Empfangers dargestellt. Urn den Einschwingvor- 
gang zu verkOrzen, wird ein spezieller Entzerrer (TDEQ = Time Domain Equalizer) vor der 
Demodulation angeordnet. Der Block TDEQ besteht hierbei aus einem Transversalfilter, 
deren Koeffizienten adaptiv eingestellt werden (siehe Bild 4) . Man spricht hier von einem 
Zeitbereichsentzerrer, obwohl hier nicht die linearen Verzerrungen des Kanals ausgeglichen 
werden. Die Aufgabe dieses Ze'rtbereichsehtzerrers ist lediglich, die Impulsantwort des Ge- 
samtsystems bestehend aus Sender , Obertragungskanal und Zeitbereichsentzemer zu ver- 
kOrzen ( Anzahl der Impulsantwortwerte mSglichst < M ), ohne zusatzliche Sperrstellen im 
Durchlaabereich zu erzeugen. Auf das Rauschverhalten hat diqser Zeitbereichsentzerrer nur 
sehr geringen EinfluB, da der Frequenzberelchsentzerrer auch ,die Frequngangverzerrungen 
des Zeitbereichsentzerrers kompensieren muB. 

In der beigefugten Anlage ( Auszug aus einer Studie, die von Siemens in Auftrag gegeben 
wurde) sind einige Verfahren zum Entwurf derartiger Zeitbereichsentzen"er angegeben. 
Hierbei handelt es sich urn Transversalfilter, die mit der Abtastfrequenz arbeiten ( bei ADSL 
mit ungefahr 2.2 MHz) , Bei einer Filterlange von etwa 20 bis 40 Koeffizienten sind daher 
ungefahr 50 bis 100 x 10 6 Multiplikationen je Sekunde durchzufOhren. ZusStzlich ergibt slch 
ein erheblicher Aufwand for die adaptive Koeffizienteneinstellung. 
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Bild 3: Blockschaltbild eines DMT EmpfSngers 



T — r 



Bild 4: Biockschattbild des Zeitbereichsentzerrers ( Block TDEQ ) 



3. In welcher Weise lost die vorliegende Erfindungsmeldung dieses Problem? 

In der vorliegenden Erfindungsmeldung wird fur den Zeitbereichsentzerrer anstelle des ad- 
aptiven Transversalfilters mit etwa 20 bis 40 Koeffizienten ein Digitalfilter mit nur wenlgen 
festen Koeffizienten vorgeschlagen. Zur Vermeidung von "echten" Multipllkatlonen werden 
Filter mit ganzzahligen oder mit einfachen Shiftoperationen darstellbare Koeffizienten aus- 
gewiahlt 

Es zeigt sich, da(i zur VerkOrzung der Systemimpulsantwort der Zeitbereichsentzerrer eine 
Nullstelle bei der Frequenz 0 haben solfte. AuBerdem soilte der Entzerrer keine Nullstelle im 
Obertragungsfrequenzbereich besitzen, urn nicht die Obertragungskapazltat zusStzIlch ein- 
zuschranken. Daher wird der Zeitbereichsentzerrer HochpaRergenschaften aufweisen. 
Filter mit der oben beschriebenen Eigenschaft lassen sich wie folgt beschreiben: 
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a) nlchtrekureive Filter 



(D 



b) rekursive Filter 



flWe*(*) 



Re* 



mit 



c 0 - ±(} - 2' L ) L ganzzahlig 



(l-co^-'r 

Wahlt man fQr m=n dann erhalt man ein Filter mit der Obertragungsfunktion 



(2) 



(3) 



Die digitalen Rlter lassen sich somit in einfacher Weise ohne MuKiplizierer realisieren. 
Im Bild 5 ret das Blockschaltbild fur ein Filter mit der Obertragungsfunktion nach (3) darge- 
stellt Hier wurde for den Koeffizienten c 0 ein posrtiver Wert angenommen. Bei Wahl emes 
negativen Koeffizienten Sndem sich lediglich die belden Vorzeichen fQr die rQckgekoppelten 
Werte Im ersten Summierer. Wahlt man fQr L=0, dann erhalt man fQr den Koeffizienten c 0 =0 
und es ergibt sich ein rein nichtrekursives System. 




insgesamt n Stufen 



i 




Bild 5: Blockschaltbild fur ein Entzerrerfilter mit der Obertragungsfunktion nach Gleichung (3) 



Bei Wahl eines mehrstufigen Filters ist es natOrlich auch moglich. die RQckkoppelkoeffizien- 
ten c 0 verschieden zu dimensionieren. Allgemein erhSlt man fur die Obertragungsfunktion 
des Zeitbereichsentzerrers den folgenden Ausdruck 
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#(*) = TT 1 * J mit =±(1-2^) 



Fur den Zeitbereichsentzerrer ergibt sich somlt ein digltales Filter, das aus einer Kettenechal- 
tung mehrerer Einzelfilter LOrdnung zusammengesetzt werden kann. Jedes Einzelfilter ISftt 
sich mit Hilfe eines Zustandsspeichers und Schiebe- sowie Additions / Subtraktlonsoperatio- 
nen realisieren. 

Urn den Realisierungsaufwand gering zu halten, ist es naturtich slnnvoll, die Anzahl der, 
Stufen und damit den Filtergrad des Digitalfilters mdglichst kleln zu halten. Nachfolgend 
werden beispielhaft die folgenden Systeme betrachtet 



Ha 



H 9 



A-05-M- 1 . 
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Die Dampfungsveriaufe der Filter H1, H2, H3 sind im Bild 6 und die der Filter H4, H5, H6 im 
Bild 7 dargestellt. Als Abtastfrequenz der Filter wird der fur den ADSL Standard spezifizierte 
Wert von 2.208 MHz gewahlt. Die Dampfungsveriaufe sind nur bis zu 025-/^ dargestellt 
FGr niedrige Frequenzen steigt die DSmpfung - abhangig von der Dimensionierung - mehr 
oder weniger stark an. Der im Empfanger nachgeschaltete Frequenzbereichsentzerrer mull 
die Abschwachung be! niedrigen Frequenzen wieder ausgleichen. Dies kann bei groBen 
Dampfungswerten zu Dynamikproblennen fOhren. Die in dieser Erfindungsmeldung beschrie- 
benen Filter sind daher dann besonders vorteilhaft einsetzbar, wenn der untere Frequenzbe- 
reich nicht fur die DatenQbertragung genutzt wird. FGr den splitterlosen digitalen Teilnehme- 
ranschlufi wir zur Zeit sowohl bei ETZI als auch bei ITU ein DMT-Verfahren unter dem 
Stichwort "ADSL-lite" diskutiert, das nur den Frequenzbereich von etwa 70 kHz bis 550 kHz 
nutzen soli. In diesem Frequenzbereich sind die hier beschriebenen Filter besonders vorteil- 
haft einsetzbar. 
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500 kHz 600 



BMJi Dampfungsverlaufe der Filter H1(z), H2(z), H3(z) 




„ 100 200 300 400 

Biid 7: Dampfungsverlaufe der Filter H4(z), H5(z), H6(z) 



kHz 



600 



Zur Beurteilung der Leistungsfahigkeit der hier beschriebenen Filter werden abschl eBend 
gVsimulationsergebnisse wiedergegeben. FOr ein ADSL Ubertragungssystem mrt *ner 
iSattleitung der Lange 3.0 km und dem Durchmesser von 0.4 mm wurden mil H.lfe «- 
£?£S2^ S/NA/erhaltnisse am Entscheidereingang for ahe *m^^~^ 
nale f Frequenzen 7 bis 250 entsprechend einem Frequenzbere.ch von etwa 30 kHz bis 1.1 
MHz) bSE Dabei wird kein au&eres Rauschen eingespeist und auch ke.ne Real.s.e- 
mng ungenauigkelten von Sender und Empfanger berucksichtigt. urn den 
hVrlichsentzerrers auf den resultierenden Fehler bzw. dem s.ch ergebendem Signatetorab- 
SSd*!S2?STiSi«"tllchan. Die Bilder 8 und 9 zeigen die Ergebnisse. die m.t den oben 
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definierten Digitalfiltem erreicht werden konnen. Zum Vergleich ist jeweils der Vertauf fQr ein 
System ohne Zeitbereichsentzerrer (schlechtester S/N-Verlauf) und ein System mit einem 
auf den benutzen Obertragungskanal im Zeitbereich optimierten Zeitbereichsentzerrer mtt 32 
Koeffizienten mtt dargestellt (bester S/N-Verlauf). Bel Einsatz elnes Zeitbereichsentzerrers 
der hier beschriebenen Art ergeben sich bereits bei Systemen 1. Ordnung erhebliche Ver- 
besserungen im Vergleich zu Systemen ohne Zeitbereichsentzerrer; Bei Systemen mit einer 
Ordnung > 1 erhait man mit Ausnahme des untersten Frequenzbereichs etwa den gleichen 
Fehler wie mit einem optimierten Zeitbereichsentzerrer mit 32 Koeffizienten. 




Bild 8: Signalstorabstand am Eingang des Entscheidere fur alle genutzten Frequenzkanele, 
ohne auBeres Rauschen und Realisierungsungenauigkeiten 




Bild 9: Signalstorabstand am Eingang des Entscheiders fQr alle genutzten Frequenzkanale, 
ohne auderes Rauschen und Realisierungsungenauigkeiten 
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1 Einleitung 



Im Rahmen dieser Studie sollten die Moglichkeiten einer im Vergleich zu ADSL (Asymmetrical* 
Digital Subscriber Line) schmalbandigeren Hochrateriubertragung erortert werden. Der vermitt- 
lungsseitige Transceiver (LT) sollte dabei als in die momentane Line-Card-Infrastrnktur der di- 
gitalen Vermittlungseinheit integriert angenommen werden. Dies resultierte in Anforderungen fur 
die mittlere Sendeleistung und bedingte eine vergleichsweise hohe, durch die Schaltungsajiordnung 
bedingte, Rauschleistungsdichte des Hintergrundrauschens. Die zugrundeliegenden "Oberlegungen 
(von Seiten des Auftraggebers) zur Kombination des 'High-Speed-Modems' mit Einrichtungen der 
Vermittlungsstelle werden hier nicht dargestellt, sind aber beispielsweise einem Konferenzpapier 
[1] zu entnehmen. 

Von dem Auftraggeber wurden Aussagen inbesondere zu drei Bandbreitestufen (32 kHz, 128 kHz, 
ca. 500 kHz) fur wichtig erachtet. Es sollten Untersuchungen'zu moglichen Ubertragungsraten 
durchgefuhrt werden. In dem vorliegenden Abschlufibericht werden einmal theoretische Ergebnisse 
auf der Grundlage von Kaualkapazitatsberechnungen, aber auch Resultate von Simuiationslaufen 
mit einer DMT-(Digcrete-Multitone)-Mehrtr^eriibertragung dargestellt. 

^Um uberhaupt zu Aussageh-gelangen^iiT^nnenT^aufiten im interessi^D^^T^^^Serm^^ 
Messungen und Modellierungen der Kabel- und der Nebensprecbubertraguiigsfanktionen durch- 
gefuhrt werden- Zusammen mit realistischen Stonumgebungen, die detailliert erlautert werden, 
bietet dies die Grundlage fur alle Berechnungen und Simulationen. 

Ein wichtiger Bestandteil dieser vorliegenden Arbeit bildef auch die Darstellung der'konkuxrieren- * 
den Ubertragungsverfahren. Es wird hierzu ebenfalls angegeben, wo diese..Verfahren Eingang in 
die Standardisierung gefunden haben, welche Herstellerfirmen sich fiir welche Verfahren stark ma- 
chen und wie sich die Schutzrechtslage gestaltet. Die Multiplexveifahr en. werden ebenso erlSutert 
und auch in kurzer Form ein Einblick in Moglichkeiten fehlerkorrigierender Codienmg gegeben. 

Da die Studie international aktuell in den Bereich der Bestrebungen zur Definition eines ADSL- 
Lite-Systems (G.lite, UAWG) gehort, wurde Uberlegungen dizu ein eigener Abschnitt gewidmet. 



Auftraggeber: Dr.-Ing. Rudi Frenzel, Siemens AG, Bereich Halbleiter, Miinchen 



Bearbeiter: Dr.-Ing. Werner Henkel und Dipl.-Ing. Thomas Kefiler, Deutsche Telekom 
AG, Technologiezentrum Darmstadt 



9. WHILE (Rz < R T ) { 

suche die Tragerposition i mit maximaiem A-fl,- 

(nur diejenigen i, fttr die gilt: Ri < Rm**) 

Ri := Ri + 1 

ARit^&Ri-l 

IF (Ri - 1) N wd ■- Nut* + 1 

} 

10, Berechnen der f% fur alle verwendeten TVager bei Vorgabe der geforderten Bitfehlerwahr- 
scheinlichkeit TV Bei quadratischen QAM-Konstellationen mit Ri > 4 ist die Leistung pro 
TVager Pi naherungsweise gegeben durch [15]: 

11. Berechne Pz = Hi,R t> o Pi u»d korrigiere die Teilleistungen entspreehend ^ 

Pt 



lm Algorithmus nach Fischer et aL wurde eigentUch die Quantisierung erst vor der Berechnung 
des QuantisierungsfeUers in Schritt 6 durchgefuhrt. Die Formulierting hier durfte sich jedoch als : j§ 
vortcilhafter herausstellen. 1*3 

' ' ' 

3.4.2 Zeitbereichsentzerrung bei DMT 

Wir hatten schon oben angesprochen, dafi eine einfache Entzerrung im DFT-Bereich durch AGC .4 
fur jeden Trager nur darm moglich ist, wenn die Kanalimpulsantwort auf das Guard- Intervall 
beschrankt bleibt. Es mufl also sichergestellt sein, dafi die Inter-Block- Interferenz moglichst gering 
ist. Eine Zeitbereichsvorentzerrung ist iiblicherweise erforderlipfcu 

Es existieren momentan drei verschiedene Ansatze zur Einstellung des Zeitbereichsentzerrers: ~ 

1. Gebrochen-rationale Kanalapproximation, \\ 

2. Optimierung des quadratischen Fehlers relativ zum Ausgang eines Ersatzsystems aus Kanal 
und Entzerrer, 

3. Optimierung der Kanalkapa2itat unter Berucksichtigung aller Storanteile (extern und Inter- 
Block-Interferenz) und des Leck-Effekts der DFT 

Die naheliegendstc Vorgehensweise ist die erstgenannte, wobei die Parameterwahl nicht immer 
von vornherein auf die Qualitat der Ergebnisse schliefien laflt. Der in der praktischen Realisierung 
verbreiteteste Abgieichalgorithmus ist sicher der zweitgenannte. Dieser fuhrt auch i.allg. auf recht 
gute Ergebnisse. Der drittgenannte Algorithmus wurde erst kiirzlich entwickelt und ist daher 
noch nicht praktisch eingesetzt worden. Simulationen zeigen jedoch, dafi er fur typische ADSL- 
Parameter, -Reichweiten und Storumgebungen im TeilnehmeranschluBnetz recht gut geeignet ist. 
Dieser Algorithmus ist in DeutscWand zum Patent (Deutsche Telekom) angerneldet. Im folgenden 
werden alle Verfahren skizziert, wobei auf das letztgenannte etwas genauer eingegangen wird. 
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1) Gebrochen-rationale Kanalapproximation: Hier wird die Kabelubcrtragurigsfunktion 
durch eine gebrochen-rationale Funktion 

H{z) = A{z)/B{z) (39) 

approximiert. Ein Entzerrer mit der Ubertragungsfunktion W(z) = B(z) reduziert die Uber- 
tragungsfunktion auf A(z), d.h. der Nenner der Kanaliibertragungsfunktion wird eliminiert. Die 
gebrochen-rationale Approximation kann relativ leicht beispielsweise mit einer Matlab-Routine 
(invfreqz der Signal-Processing Toolbox) durchgefuhrt werden und fuhrt auch, wena man die 
Ordnung der ZaJiler- und Nennerpolynome etwas variiert, zu sinnvollen Ergebnissen. Hierdurch 
wird jedoch nur die Impulsantwort verkiirzt und die Storumgebung iiberhaupt nicht beriicksich- 
tigt Audi ergeben sich teils bessere Ergebnisse, wenn man den Zahlergrad grofler ale das Guard- 
Intervall wahlt, obwohl dies eigentlich die Design-Regel sein oolite, denn das Zahlerpolynom kenn- 
zeichnet ja die Restimpulsantwort. Dies hat wohl numerische Ursachen und hangt sicher von den 
verwendeten Approximationsalgorithmen ab. 

Das Verfahren der gebrochen-rationalen Approximation wurde in dieser Studie nur einmal alter- 
nativ bei der sehr schmalbandigen Anwendung mit einer Bandbreite von nur 32 kHz eingesetzt, da 
sich hie rder-neue Algorithmus 3, wie er moment an realisiert ist r als-weniger geeignet heraus ges tellt^= 
hat. 

2) Ersatzsystem-Methode: Alle Vorschlage, die hierzu existieren, gehen wohl im Prinzip auf 
eine Veroffentlichung von Falconer und Magee, 1973 [16], zuruck. VerofFentlichungen neueren Da- 
tums sind beispielsweise die von Bladel, Moeneclaey [17] und Faxhang-Boroujeny [18], Einleitend 
ist der Ansatz jedoch auch in einer Arbeit von Al-Dhahir, Cioffi [19] kurz mit beschrieben. Auf 
diese Arbeit wird bei der Erlauterung des dritten Abgleichalgorithmus nochmals hingewiesen. 

Bild 33 zeigt die prinzipielle Anordnung auf die sich alle diese VorschlSge beziehen. Man stellt 
sich hier ein Ersatzsystem zu der Kettenschaltung von Kanal und Entzerrer vor. Der quadrati- 
sche Fehler zwischen den Ausgangssignalen von Kanal/Ent^errer und Ersatzsystem ist zu mini- 
mieren. Dies geschieht durch wechselseitiges Anpassen der Koeffizienten des Entzerrers und des 
Ersatzsystem^. Die Herleitung des grundsatzlichen Verfahrens soil hier kurz dargestellt werden. 
Es existieren mannigfache Simplifikationen, die auch beispielsweise in einem Patent von Amati 
aufgefuhrt werden. Auf diese Varianten werden wir hier nicht detailliert eingehen. Oblicherweise 
werden dort Teilschritte zur Optimierung durch Abschneidoperationen ersetzt, d.h, beispielsweise 
kann ein zulassiges Ersatzsystem einfach durch Langenbegrenzung der Gesamtimpulsantwort ge- 
wonnen werden. Ahnliches gilt fur den Entzerrer. Die Optimierung beschrankt sich dann nur noch 
auf eine Komponente: Entzerrer oder Ersatzsystem. Man konnte sogar den gaiizen Algorithmus 
durch solche Abschneidoperationen realisieren, jedoch waren eigene Untersuchungen liierzu wenig 
erfolgversprechend. 

Zur genaueren Beschreibung des Algorithmus gehen wir von einer vektoriellen Beschreibung aus. 
h = [hoj. . . ,/iA*-i] r bezeichne die urspningliche Kanalimpulsantwort ohne Entzerrung, w s= 
[v; 0 ,. % w nw -\] T sei der Entzerrer und b = [6 0 , . . . , 6 n/# _,] r sei die gewunschte (kurze) Impuls- 
antwort des iSrsatzsystems, das die Kettenschaltung von h und w reprasentiert. v sei ein Vektor 
aus n w empfangenen Abtastwerten und x 6 sei ein Vektor aus n 6 Abtastwerten des Sendesignals, 
allerdings urn A verzogert. Unser 2iel ist es nun, w, b und A so zu bestimmen, dafi der mittlere 
quadratische Fehler 7,w]$chen w r • v und b r • x& minimiert wird. 



33 



Kan a) 




VerzOgemng Ereatzsystem 
Bild 33: Ersatzsystem von Kanalimpulsantwort und Zeitbereichseatzoref 

Der mittlere quadratische Fehlcr kann wie folgt geschrieben werden: 

E{e'} » £{(w T .v-b r .x*) 2 } 

= w r R w w + b T R«b - 2w T R.xb 



(40) 



mit 



Raw = £?{Xfc xf} 



(fife 

(42) - 

(43) : 



Rvx und Rsa hangen beide von der Verzogerung Aab. 

Fur eine vorgegebene Verzogerung A kann die OptimaleinateUung der Entzerrerkoeffizienten 
dureh Nullsetzen des Gradienten des mittleren quadratischen Fehlers gefunden werden. Anwenden 
von dB{e 2 }/dvf = 0 fuhrt auf . 

vropt = R^IUb (44) 

Dies stellt nun das optimale w in Abhangigkeit von b dar. §etzt man dies wiederum in (40) ein, 
so erhalt man ■ • ^ , /.,-\ 

f?{e 2 } = b T .(iUx-R2L(iC 1 ) r R-) ta = br ' c - b W 

Dieser Ausdruck ist zu minimieren. Dabei ist jedoch die triviale und unsinnige Losung b = O r j 
vermeiden. Um dies sicherzustellen, wird als Nebenbedingung die Normierung b x b = 1 verwendet. 
Ein Ansatz mit Lagrange-Multiplier liefert dann 



E{e 2 } = b r .C-b + A(l-b r -b) 



dh 



= 2C • b + 2Ab = 0 
C ■ b = Ab . 



(46) 

(47) 
(48) 



Es ist somit ein Eigenwertproblem zu losen. Setzt man die LSsung in Gl. 46 ein, so erhalt man 

E{c 2 } = A < b T • b = A (49) 

Die resultierende quadratische Abweichung ist somit gleich dem Eigenwert. Damit ist die gesuchte 
Losung gleich dem Eigenvektor, der zum kleinsten Eigenwert gehort. 
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Damit wollen wir die prinzipieUe Darstellung dieses Algorithmus abschliefien. Die Beschreibung 
aller Varianten des Algorithmus wurde den Rahmen dieser Studie sprengen, 

Es ist anzunehmen, dafi alle momentan verfugbaren ADSL-Modems eine Variante des Ersatzsystem- 
Abgleichalgorithnius verwenden. Beispielsweise diirfte die oben erwahnte Veroffentlichung von Bla- 
del, Moenecjaey [17] die Grundlage der Implementationen von Alcatel bilden, da die Arbeiten zu 
dieser Veroffentlichung von Alcatel beauftragt wurden. 

3) Optimierung der Kanalkapazitat: Bei diesem Algorithmus wird die Kanalkapazitat 

fi=L^=E^ 1 + ^ - (50) 

maximiert. Die Maximierung der Kanalkapazit&t entspricht naherungsweise der Maximierung des 
geometrischeu Mittels aus den Storabstanden (1 im Argument des Logarithmus vernachlassigt). 
Es werden alle Rauscheinflusse, sowohl externe als auch durch Inter-Block-Interferenz bedingte, 
zusammengefafit, und es wird auch berucksichtigt, dafi die empfangsseitige DFT eine Rechteck- 
fensterung darstellt, die zu einem Leckeffekt fuhrt. Sehr offensichtlich wird der Leckeffekt in der 
=^ahe von DC. Dort sollte idealerweise der "Stdrabs Land gegen unendlich gehen, wenri man' von der 
NEXT-Funktion ausgeht, die dort eine Nullstelle hat. Auch die Ubertragungsfunktion ist um DC 
am gunstigsten. Tatsachlich wird jedoch der Storabstand in der Gegen d um DC wieder schlechter, 
da Rauschanteile infolge des Leckeffekts dorthin gefaltet werden (siehe die SNRrDarstellungen in 
Abschnitt 10). Der Leckeffekt ist lediglich fur das Rauschen von Bedeutung, da dieses im Gegensatz 
zum Signal keinen zyklischen Prafix aufweist. 

Der Algorithmus findet eine Losung, die ublicherweise nicht streng die Impulsantwort auf die Lange 
des Guard-Intervalls (+1) beschrankt. Er sucht stattdessen einen Kompromifl zwischen externem 
Gerausch und Inter-Block-Interferenz, wobei der Leckeffekt eine wesentliche Rolle spielt. Ohne 
den Leckeffekt ware auch lediglich die Impulsantwortlange von Bedeutung, denn eine beliebige 
Obertragungsfunktion wurde das SNR fur die einzelnen Trager nicht verandern, da ja sowohl 
Signal als auch Rauschen mit dem gleichen Faktor multiplizi^rt wurden. 

Die Uberlegung, die Kanalkapazitat zu optimieren, geht zuriick auf eine Veroffentlichung von Al- 
Dhahir, Cioffi [19], wo zwar dieses Optimierungskriterium genannt wird, bei der Dnrchfuhrung 
jedoch wird nur das geometrische Mittei der Kanalenergie maximiert. Zusatzlich wird wieder auf 
den Ersatzsystem-Algorithmus zuruckgegriffen. Die zentrale Erkenntnis, die es nun erst gestattet, 
die Kanalkapazitat selbst wirklich zu optimieren, ist das Einbeziehen des Leckeffekts der DFT, 
welcher zuvorbi slang nicht beriicksichtigt wurde. 

Der Einfachheit halber zitieren wir nun zur Beschreibung des Algorithmus eine kurzlich erstellte 
Veroffentlichung in englischer Sprache (der deutsche Patenttext ist durch unabgestimmte Ande- 
rungen der Patentabteilung etwas unverstandlich). 

In the sequel, the steps of the iterative algorithm arc described. There, N^tr denotes the maximum 
possible number of used carriers. For baseband applications like DMT, one has to consider conju- 
gacy constraints to ensure that the time-domain signal is real. This means that the upper half of 
the DFT frame is dependent on the lower half. The DFT block length is then Noft = 2 • 7V carr . 
The algorithm, however, requires an ovcrsampling in frequency domain by a factor of two, which 
leads to a DFT block length of Ndft* = 4 • A^ rr . 
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The steps of the adaptation algorithm 1 : 

1. Determine the DFT transform fsf = of the impulse response corresponding to the noise,' 
taking into account all the spectral properties of the external noise. is the discrete noise 
power spectral density. 

2. Let Hi be proportional to the discrete channel response. The quotient Hi/M{ is normalized, 
such that -"d* rW*| a /]C* W'l 3 represents the signal-to-noise ratio- SNRin at the input of the 
equalizer when using all possible carriers with the same transmit power: 



Hi : — Hi 



N SNRi -™ (S1) 



3. Product S; of the channel response %- times the equalizer response W{ in the DFT domain: 

Si^Tii Wi * (52) .j 

4- Determine the corresponding impulse response by applying the IFFT: 

3 = 1FFT(S) — " (53) ^ 

5. Determine the mean signal power density fraction \Nf\ 2 of the part of the impulse response ^jj 
s that exceeds the guard interval and influences the next frame. 

6. Compute the product of Mi times the equalizer frequency response Wi: "*f 

N?=Mi-Wi (54) ~ 

7. Add the noise power density components from steps 5 and 6: 

|/V;-| 2 = IW/I 2 (55) 

i 

8. Compute the autocorrelation function n t - that corresponds to the total noise power spectral 
density \Ni\ 2 : 

(n,.)=IPPT( (\N<\ 2 ) ) (56) 

9. Introduce the leakage effect of the DFT by multiplying n; with a triangular function 

A; = m - ti , 2 = 0,..., 4iVcarr - 1 (57) 



According to [20], the influence of a rectangular time-domain windowing, which is due to 
the DFT operation at the receiver, corresponds to the multiplication of the time-domain 
equivalent (autocorrelation function) of the noise power spectral density times a triangular 
function. In order to describe the triangular function approximately with the discrete time- 
domain vector, an oversampling in frequency domain by (at least) a factor of 2 is chosen. 
This means that the algorithm relies on a frequency scale with at least the half of the DMT 
carrier spacing. 



denotes the vector with components x f . 
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10. Traasform of A into the DFT domain: 

A = FFT(A) (58) 

11. Determine the signal-to-noise ratio SNRf for every used carrier location. No frequencies are 
considered that just resulted from the oversampling, only the ones that are actually used 
and are at multiples of the original carrier spacing. 

SNR* = S</A ; , Si « \Sif , t 6 U (59) 

U denotes the set of actually used carrier locations. 

12. Summation of all channel capacities at all carrier locations: 

-R = £log 2 (l+SNR<) (60) 

R is the parameter to be optimized, which can be carried out by any multidimensional optimization 
algorithm, like AMOEBA (DownfifH^Simplex Method) in [21], it^afively modifying tHe time- ' 
domain equalizer coefficients. For each iteration, all the steps 3 to 12 have to be repeated* The 
change of R between two successive iterations may serve as a termination criterion. This means, 
the actual value of R is stored for the next iteration to allow for a comparison, afterwards. If the 
change of R falls below a certain threshold, the iteration is terminated. 
Instead of R 7 the geometric mean of all the SNRf, t € U, can be maximized. 

Both noise components (external noise and inter-block interference) are combined in Step 7 and 
the leakage is treated in Step 9. 

Step 5 needs to be explained in some more detail. When receiving a DMT signal, the detection 
frame will be positioned in time such that the inter-block interference is minimized. The optimum 
position is a tradeoff, where precursors resulting from components of the preceding frame as well 
as postcursors resulting from components of the following frame have to be considered as a distur- 
bance. Pour such exemplary impulse responses resulting from components of neighboring frames 
next to the detection frame are depicted (for illustration purposes in continuous form) in Fig. 34, 
The rectangular dots symbolize further impulse responses resulting from the following components 
on both sides of the detection frame. For every position of the impulse responses, the portion in- 
side the detection frame has to be extracted by a rectangular window. The number of impulse 
responses that have to be considered is dependent on the length of the impulse response. Apply- 
ing an FFT and afterwards computing the squared amplitudes of every DFT-domain component 
(periodogram) yields the corresponding noise power spectral density (PSD). All such noise PSD 
contributions from the neighboring frames have to be summed up to obtain a mean noise PSD of 
the inter-block interference. Like for signal-independent external noise, the rectangular windowing 
of the DFT has to be considered. This means the leakage effect comes into play here, too. For 
illustration, both noise types are summed up in Step 7, before the leakage is incorporated in Step 
9. Hereto, the detection frame has to be taken twice as long, like for external noise. However, the 
computationally more efficient solution would be to use a non-extended detection frame, whereby 
the leakage is automatically taken into account. The less computationally efficient description has 
been prefered here, because it unveils the similarity between both noise components more clearly. 
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Let the length of the guard interval be n 0 . The position of the detection frame can be estimated 
by first finding the position of a frame of length no + I with the maximum energy content when 
sliding over the entire impulse response. The detection frame starts right after it. 

In the sequel, the procedure is again described in a more formal way. Let us first suppose that 
the time axis zero point of the total impulse response 5,- is located at the 'center of energy \ Such 
a time-axis shift has no influence on the results. Now, let the impulse response have non-zero 
components for i = -n v , . . . , 0, . . . , n m with n v being the number of precursor and n„ being the 
number of postcursor samples. The mean PSD portion | Nf | 2 resulting from the part of the impulse 
response that influences the neighboring frames can be computed as follows: 

Precursors: Va u = •0,...,n u - a (possible positions of the impulse response, resulting from 
different samples; a, see, Fig. 34): 



S 4N^n—l-i = 5 — a-i— er„ , t = 0, . . . , Tl v — Q — <T. 
£ <LN CB TT-l-i = 0 i i = »»v — a — C v + 1, 



V 



Post cursors ?= .0, . " n — nxr-R t.— 2: ^ • - 

^n-na + a- 2 — cr n * 

As already described, we chose the frame length of the detection frame to be 4N cm samples - 
long, because of similar treatment of both noise components.. The real detection frame at a DMT 
receiver has, of course, only 2N CAtr samples. 

The mean PSD of the inter-block interference is now obtained by summing up all precursor and 
postcursor components. 

with - FFT(s<°")) , = FFTfal**)) . .;' 'jjj' 

Bei der Anwendung des Algorithmus fur die sehr geringe Bandbreite von 32 kHz ergaben sich bei 
der Adaption Stabilitatsprobleme. Dies liegt jedoch sicher nicht in dem Prinzip des Algorithmus 
begriindet, sondera wahrscheinlich in der programmtechnischen Realisierung. Bei den grofleren 
Bandbreiten wurden die Parameter (Abtastrate, Tragerabstand) von ADSL zugrundegelegt. Dort 
gab es keinerlei Probleme mit dem Entzerreralgorithmus 3. Urn zu klaren, warum der neue Entzer- 
reralgorithmus bei einer Bandbreite von nur 32 kHz wenig geeignete Entzerrereinstellungen lieferte, 
ware eine langere Untersuchung notig gewesen, die im Rahmen der Studie nicht zu leisten war. 
Mit einer geringfugigen Modification im Algorithmus 3 wurde dennocli eine Entzerrereinstellung 
gewonnen. Es zeigte sich jedoch auch hier ein instabiles Verhalten wahrend der Adaption. 

Dies schliefit die Behandlung der beiden zentralen Algorithmen - Bitallokation und Zeitbereichs- 
entzerrereinstellung - eines DMT-Systems ab. Wir hatten schon oben angemerkt, dafi eine Proble- 
matik eines DMT-Systems das relativ hohe Spitzen-Mittelwert-Verhaltnis ist. Dieser Problematik 
ist das folgende Kapitei gewidmet. 
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Bild 34: Inter-Block-Interferenz von Nachbarrahmen 



3,4-3 Reduktion des Spitzen-Mittlerwert-Verhaltmsses bei DMT 

Das hohe Spitzen-Mittlerwert-Verhaltnjs (Peak-to- Average Ratio, PARsCrest-Faktor 2 ) resultiert 
aus der nahezu Gaufl-formigen Amphtuden-H&u£gkeitsverteilung eines DMT-Zeit signals. Dies 
steilt besondere Anforderungen an die Dynamik der analogen Schaltungsteile. In einem spatereix 
Kapitel wird aus dem Spitzenwert, den die Analogschaltung verarbeiten kann, die mittlere Sen- 
deleistung des High-Speed-Modem berechnet. Dabei liegt die durch das Abschneiden (Clipping) 
hervorgerufene Storleistung und die daraus resultierende Bitfehlerwahrscheinlichkeit zugrunde. Im 
Falle von Echokompensation sollen die Auswirkungen wegen nichtlinearer Echos noch kritischer 
sein. Dazu liegen jedoch bei uns bislang keine Erfabirungen vor. 

Es existieren bislang funf prinzipiell verschiedene Vorschlage zur Reduzierung der Spitzenwer- 
te, wobei insbesondere einer hervorzuheben ist, da er vollstandig im Zeitbereich arbeitet und 
damit einen relativ geringen Aufwand bedeutet. Algorithmen, die mehrere FFTs oder gar DFT- 
Transformationen erfordern, sind praktisch fur eine Anwendung bei xDSL nicht von Interesse. Al- 
ternativ ware vielleicht auch ein Verfahren denkbax, das nur im DFT-Bereich operiert. Allerdings 
existiert hierzu nur ein Vorschlag (Patent Deutsche Telekom), der aus Grunden des Speicherbe- 
darfs sich nur fur 4-^PSK auf alien Tragern und relativ geringen Tragerzahlen (< 64) eignet. 

Es wird nun kurz eine Ubersicht iiber die verschiedenen Verfahren gegeben. Das momentan gunstig- 
ste, bercits angesprochene Verfahren, das auch bei ITU Study Group 15 als Beitrag [22] eingebracht 
wurde, wird ctwas ausfuhrlicher dargestellt. 

Iterative Spitzenwerteliminierung durch Projektion auf redundante Unterraume [23] 
Dieses Verfahren verwendet (aquidistante) nicht belegte Trager als Redundanzstellen, mit denen 
die abgeschnittenen Signalspitzen eliminiert werden. Die Redundanzstellen spannen einen Unter- 
raum auf, auf den der Signalverlauf aus abgeschnittenen Signalanteilen projiziert wird. Es wird der 
entsprechend negative Wert den redundanten Tragerbelegungen hinzuaddiert. Nach einer Quanti- 
sierung der Tragerbelegung auf die zulassigen Werte erfolgt eine Iteration, die allerdings jedesmal 
die Berechnung des Zcitbereichssignals enlhalt. 
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Notification of invention 



Circuit arrangement for a time domain equalizer for discrete multitone 
modulation 

1. What technical problem is the invention intended to solve? 

Discrete multitone modulation (DMT) can be used to transmit data over channels 
which introduce linear distortion. In this case, advantages may arise over single- 
carrier methods in terms of flexibility when matching the data rate or the 
transmission spectrum to the transmission channel or the noise environment. The 
method Is used for. digital broadcasting (DAB - Digital Audio Broadcast) under the 
name OFDM (Orthogonal Frequency Division Multiplex). This method has likewise 
been standardized for the asymmetric digital subscriber line (ADSL) using normal 
telephone lines. 

In the modulation method, the transmission signal is composed of a large number of 
sinusoidal signals, with each individual sinusoidal signal having both its amplitude 
and its phase modulated. What is obtained is therefore a large number of 
quadrature amplitude modulated signals. 

This is implemented by using inverse FFT (Fast Fourier Transformation) in the 
transmitter and normal FFT, for which there are efficient signal processing 
algorithms, in the receiver. « 

H 

Figure 1 shows the basic block diagram for the DMT transmission system. 




Noise 



HL DT CE 
Dr. H. Schenk 



Drawn up on: 
Printed on: « 



ZEITE2_1.doc 



-2- 



Fioure 1 : Block diagram of the DMT transmission system. 

The data to be transmitted are first combined into blocks. Depending on hierarchy, a 
certain number of bits to be transmitted is then respectively allocated to a complex 
number. N such complex data vectors are then subsequently transformed Into the 
time domain using inverse discrete Fourier transformation (IDFT). The values 
obtained after transformation are immediately the samples of the transmission 
signal. The complex data vectors may thus be regarded as amplitudes for cosine 
(real part) and sine (imaginary part) oscillations which are to be transmitted within a 
block. The frequencies are distributed equidistantly in this case and are 

fl = l .l 2, .....N/2 

where T is the duration of a block and N is the number of samples in a DFT block. If 
a power of two is chosen for N, the inverse FFT (Fast Fourier Transform) can be 
applied for the transformation into the time domain, in which case the 
implementation complexity is severely reduced. 

Before being forwarded to the D/A converter, the last M samples of a block (M<N) 
are appended to the start of the block again (cyclic prefix). This allows a periodic 
signal to be simulated in the receiver when the transient process caused by the 
transmission channel has died down after M samples, in line with a time of T-— . 

The complexity of equalization can thus be severely reduced in the receiver, since 
after the demodulation in the receiver it is only necessary to multiply by the inverse 
channel frequency response in order to eliminate the linear distortions In the 
transmission channel (one complex or four real multiplication operations per 

frequency line). . 
In the case of data transmission over copper lines, relatively long transient 
processes arise. On the other hand, it is desirable to keep the cyclic prefix as small 
as possible in relation to the block length (M«N) in order to obtain the smallest 
possible reduction in the transmission capacity. By way of example, in the ADSL 
standard a block length of N-512 has been stipulated for the data transmission from 
the exchange to the subscriber and a value of M=32 has been stipulated for the 
cyclic prefix. Since the transient process has not yet died down after 32 samples, 
additional errors arise in the receiver which cannot be compensated for by the 
frequency domain equalizer. Figure 2 shows the S/N ratio in dB, obtained using a 
simulation, after the demodulation and ideal frequency domain equalization directly 
upstream of the decision-maker for all frequency channels used for an ADSL 
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system. Of the N=256 frequency lines, lines 7 to 250 are used for data transmission. 
The line interval is 4.3125 kHz, which means that the range from approximately 
30 kHz to 1.1 MHz is obtained for the transmission frequency band. The 
transmission line chosen was a two-wire line with a length of 3 km and a diameter of 
0.4 mm (parameters based on ETSI DTR/TM-06001 ). 

Ignoring external and internal interference (channel noise, crosstalk, implementation 
inaccuracies for the A/D and D/A converter and for the transmission amplifier, 
quantization noise), values in the range from approximately -21 to -38 dB are 
obtained for the S/N ratio at the decision-maker input. 

The error caused by the transient component not dying down after the demodulation 
can be reduced or eliminated using special signal processing measures. The 
present notification of invention describes a method and a circuit arrangement for 
eliminating such interference. The implementation complexity measured as the 
number of multiplication operations per second for the correction and for the 
adaptive adjustment is much lower in the case of this method than in the case of the 
methods used at present. 




Fioure 2: Signal-to-noise ratio at the input of the decision-maker for all frequency 
channels used, without a time domain equalizer, with an ideal 
frequency domain equalizer, line length 3 km (0.4 mm diameter) 

2. How has this problem been solved to date? 

Figure 3 shows the block diagram of a DMT receiver. To shorten the transient 
process, a special equalizer (TDEQ ■ Time Domain Equalizer) is arranged upstream 
of the demodulation section. In this context, the block TDEQ comprises a transverse 
filter whose coefficients are adjusted adaptively (see figure 4). This is referred to as 
a time domain equalizer, although the linear distortions in the channel are not 
equalized in this case. The task of this time domain equalizer is merely to shorten 
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the impulse response of the overall system comprising transmitter, transmission 
channel and time domain equalizer (number of impulse response values as far as 
possible < M), without producing additional stop points in the passband. This time 
domain equalizer has only very little influence on the noise response, since the 
frequency range equalizer also needs to compensate for the frequency response 
distortions in the time domain equalizer. 

The attached annex (extract from a study commissioned by Siemens) indicates a 
few methods for designing such time domain equalizers. These are transverse fitters 
which operate at the sampling frequency (in the case of ADSL at approximately 
2.2 MHz), For a filter length of approximately 20 to 40 coefficients, it is therefore 
necessary to perform approximately 50 to 100 x 10 6 multiplication operations per 
second. In addition, considerable complexity is obtained for the adaptive coefficient 
adjustment 
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Figure 3: Block diagram of a DMT receiver 
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Figure 4: Block diagram of the time domain equalizer (block TDEQ) 



3. How does the present notification of invention solve this problem? 



Instead of the adaptive transverse filter with approximately 20 to 40 coefficients, the 
present notification of invention proposes a digital filter with Just a few fixed 
coefficients for the time domain equalizer. To avoid "genuine" multiplication 
operations, filters are selected which have integer coefficients or coefficients which 
can be represented using simple shift operations. 

It is evident that to shorten the system impulse response the time domain equalizer 
should have a zero point at the frequency 0. In addition, the equalizer should not 
have a zero point in the transmission frequency range, so that the transmission 
capacity is not additionally restricted. For. this reason, the time domain equalizer will 
have high-pass filter properties. Filters with the property described above can be 
described as follows: 

a) nonrecursive filters: 

b) recursive filters 

H_ = - fo Rck(Z) , with c 0 -lk-2~ L )L integer (2) 

I I 

■ t t 

If m=n is chosen, then a filter is obtained with the transfer function 



1 - Co • z" : 



(3) 



The digital filters can thus easily be produced without a multiplier, 
Figure 5 shows the block diagram for a filter with the transfer function shown in (3). 
In this case, a positive value has been assumed for the coefficient co. When a 
negative coefficient is chosen, only the two arithmetic signs for the feedback values 
in the first summator change. If L=0 is chosen, then Co=0 is obtained for the 
coefficient, and a purely nonrecursive system is obtained. 




Figure 5: Block diagram for an equalizer filter with the transfer function shown in 
equation (3) 

If a multistage filter is chosen, it is naturally also possible to provide the feedback 
coefficients Co with different dimensions- Generally, the following expression is 
obtained for the time domain equalizer's transfer function: 

Hz) = n f 1 " Z ~*"il ^th c v = ±CL - 2T*-> 

i-t V 1 - c v * 2 J 

For the time domain equalizer, a digital filter is thus obtained which can be made up 
of a concatenated circuit comprising a plurality of 1st-order single filters. Each single 
filter can be implemented using a state memory and shift and also 
addition/subtraction operations. » 

To keep down the implementation complexity, it naturally makes sense to keep the 
number of stages and hence the filter grade of the digital filter as small as possible. 
The text below considers the following systems by way of example: 
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U- 0:5-2-' J 

The damping profiles of the filters H1, H2, H3 are shown in figure 6, arid those of the 
filters H4. H5, H6 are shown in figure 7. The sampling frequency chosen for the 
filters is the value of 2.208 MHz which is specified for the ADSL standard. The 
damping profiles are shown only up to 0.25 • f A . For low frequencies, the damping 
rises - depending on the design - to a greater or lesser extent. The frequency range 
equalizer connected downstream in the receiver needs to compensate for the 
attenuation at low frequencies again. This can result in problems with dynamics at 
high damping values. The filters described in this notification of invention can 
therefore be used to particular advantage when the lower frequency range is not 
used for data transmission. For the splitterless digital subscriber line, both ETSI and 
ITU are currently discussing a DMT method under the headword "ADSL-lite", which 
is intended to use only the frequency range from approximately 70 kHz to 550 kHz. 
In this frequency range, the filters described here can be used to particular 
advantage. 




Figure 6: Damping profiles for the filters H1(z). H2(z), H3(z) 




Figure 7: Damping profiles for the filters H4(z), H5(z), H6(z) 

To assess the performance of the filters described here, a few simulation results are 
finally shown. For an ADSL transmission system with a two-wire line of length 3.0 km 
and diameter 0.4 mm, a simulation has been used to determine the S/N ratios at the 
decision-maker input for all frequency channels used (frequencies 7 to 250 In line 
with a frequency range from approximately 30 kHz to 1.1 MHz). In this case, no 
external noise is supplied and also no implementation inaccuracies for the 
transmitter and receiver are taken into account in order to provide better clarification 
of the influence of the time domain equalizer on the resultant error or on the 
resultant signal-to-noise ratio. Figures 8 and 9 show the results which can be 
achieved with the digital filters defined above. For the purposes of comparison, the 
profile for a system without a time domain equalizer (worst S/N profile) and a system 
with a 32-coefficient time domain equalizer optimized in the time domain to the 
transmission channel used (best S/N profile) is respectively shown too. When using 
a time domain equalizer of the type described here, significant improvements over 
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systems without a time domain equalizer are obtained even with let-order systems. 
In the case of systems whose order is > 1 , approximately the same error as with a 
32-coeffictent optimized time domain equalizer is obtained, with the exception of the 
lowest frequency range. 




XH2 



Figure 8: SjgnaMo-noise ratio at the input of the decision-maker for all frequency 
channels used, without external noise and implementation inaccuracies 
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Figure 9: Signal-to-noise ratio at the input of the decision-maker for all frequency 
channels used, without external noise and implementation inaccuracies 



1 Introduction 

This study will discuss the possibilities of more narrowband transmission at high 
rates in comparison with ADSL (Asymmetrical Digital Subscriber Line). In this case, 
the exchange-end transceiver (LT) should be assumed to be integrated Into the 
current line card infrastructure of the digital switching unit. This resulted in demands 
for the mean transmission power and was the cause of comparatively high noise 
power density for the background noise as the result of the circuit arrangement. The 
underlying considerations (on the part of the commissioning party) for combining the 
-high speed modem" with devices at the exchange are not Illustrated here, but can 
be found in a conference paper [1 ] for example. 

The customer has deemed statements particularly relating to three bandwidth stages 
(32 kHz 128 kHz, approximately 500 kHz) to be important. Examinations regarding 
possible transmission rates should be carried out For one thing, this final report 
illustrates theoretical results on the basis of channel capacity calculations, and also 
results from simulation passes using DMT (Discrete Multltone) multicarner 
transmission. 

To be able to obtain any statements at all, it was necessary to perform 
measurements and modeling of the cable and crosstalk transmission functions in the 
frequency range of interest. Together with realistic interference environments, which 
are explained in detail, this provides the basis for all calculations and simulations. 

An important part of this work is also the presentation of the competing transmission 
methods This is done by likewise indicating where these methods have gained entry 
to the standardization, which firms of manufacturers are putting their efforts Into 
which methods and what the situation is with industrial rights. The multiplex methods 
are likewise explained and there is also a brief insight into options for error^orrective 
coding. 

Since the study is currently internationally part of the area relating to the efforts to 
define an ADSL-lite system (G.lite, UAWG). considerations in this regard have a 
separate section devoted to them. 



Commissioned by: Dr.-ing. Rudi Frenze l, Siemens AG. Semiconductors Division, Munich 
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9. WHILE (Ri < Rt) { 

look for carrier position i with maximum ARi 

(only those i for which the following applies: R t < ) 

IF (&i =l)Nused : ^Nwxxt +I 

} 

10. Calculate Pj for all carriers used when prescribing the demanded bit error 
probability Pb. In the case of quadratic QAM configurations with Rj St 4, the 
power per carrier Pi is approximately given by [1 5]: 



2 N ± (2* 1 - 1) 



erf c" 



^v? 7 - ij 



(38) 



11. Calculate p £ = >0 p a and coiTect the power components according to 

In the algorithm according to Fischer et al., the quantization was actually not 
performed until prior to the calculation of the quantization error in step 6. In this 
case, however, the wording ought to be found to be more advantageous. 

3,4.2 Time domain equalization in DMT 

We had already mentioned above that simple equalization in the DFT range using 
AGC for each carrier is possible only if the channel impulse response remains 
limited to the guard interval. It is thus necessary to ensure that the interblock 
interference is as low as possible. Time domain preequalization is normally required. 

There are currently three different approaches to adjusting the time domain 
equalizer: 

1 . fractional-rational channel approximation, 

2. optimization of the square error relative to the output of a substitute system 
comprising a channel and an equalizer, 



3, optimization of the channel capacity taking into account all interfering 
components (externally and inter-block interference) and the leakage effect of 
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the DFT 



The most obvious practice is .he flrs.-men.lor.ed. where the >»™""""*™ 
does not alweys allow the quality of the results to be Inferred from the outset The 
Inment aigoritnm which is most widespread in the practical Implement ,s 
rtZ *e second one mentioned. This aiso generally leads to very jod**. 
The third algorithm mentioned was developed only recently and has therefore not 
y* been used In practice. However, simulations show the. it is very well suited to 
CcTaDSL parameters, ranges and interference environments in the subscriber 
Z netork. This algorithm is the subject of a patent application (Deutsche 
Tetekom) in Germany. The text below outlines all methods, the test one mentioned 
being discussed in somewhat more detail. 

,) Fracttonal-raacnal channel approximation: in this case, the cable .ransfer 
function is approximated by a fractional-rational function 

H(z) = AteVftz). < 38 > 
An equaiizer with the transfer fundon W(z) = B(z) reduces the transfer funcbon to 
Ate) I e the denominator of the channel transfer function is eliminated. The 
tSionalJttonal approximation can be performed relatfvely easily, for example 
Xatatb routine (Invfreqz from the signal processing toolbox), and produces 
usTm results even if the order of the enumerator and denominator polynomials « 
£2 "m-hat This on,y shortens the impulse response. 

.eke into account the Interference environment at all. There are also partly better 
es u robZed r the numerator degree is chosen to be greater man me guard 
nerval although .his should ectually be .he design rule, because .he numerator 
^mia, characterizes .he residua, impulse response, of course Th« probacy 
n^sTumerical causes and is certainly dependent on the a P prox,ma..on algonthme 

used. 

Th. mamod of fractional-rational approximation has been used only once In this 
Tdy af an a^ate for me very narrowband applied wKh • I-"**-* - ^ 
32 kHz, since the new algorithm 3, as is current* implemented, has been found to 
be less suitable in this case. 

2) Substitute-system method: All proposals which exist in this regard can 
1 robabt oe attributed in principle to the publication by Falconer and Magee 9173 
M6, Publications of more recent date are those by Bladel, Moeneciaey t 17land 
Famang^oroujeny [18] for example. By way of introduction, however, the approach 
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is also briefly described In a piece of work by Al-Dhahlr, Cioffi [19]. Reference will 
again be made to this work In the explanation of the third alignment algorithm. 

Figure 33 shows the basic arrangement to which all of these proposals relate. What 
is seen here is a substitute system for the concatenated circuit comprising channel 
and equalizer. The square error between the output signals from the 
channel/equalizer and the substitute system can be minimized. This is done by 
reciprocally aligning the coefficients of the equalizer and of the substitute system. 
The derivation of the basic method will be illustrated briefly in this case. There are 
multifarious simplifications which are also listed in a patent by Amati. for example. 
We will not discuss these variants in detail here. Normally, substeps for optimization 
using clipping operations are substituted there, i.e. an admissible substitute system 
can be obtained simply by limiting the length of the total impulse response, for 
example The situation is similar for the equalizer. The optimization Is then limited 
only to one component: equalizer or substitute system. It would even be possible to 
implement the entire algorithm using such clipping operations, but separate 
examinations in this respect were less promising. 

To describe the algorithm in more detail, we will use a vectorial description. Let 
h - rh„ hu.i] T denote the original channel impulse response without equalization, 

w - wl wlf be the equalizer and b = [b 0 b^f be the desired (short) 

impulse response of the substitute system, which represents the concatenated 
circuit of h and w. Let v be a vector comprising n w received samples, andoo, be a 
vector comprising n b samples of the transmission signal: but delayed by A. Our aim 
is now to determine w, b and A such that the mean square error .s minimized 
between w T • v and b T • Xb. 



Channel 



Equalizer 




Delay 



Substitue 
system 



Figure 



33: Substitute system for channel impulse response and time domain 

equalizer 



The mean square error can be written as 
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B{e*} = ^{(w T .v-b r .x 6 ) 2 } (40) 



with 

Rw = E{v • v r } . (41? 

Rvx = E{v • xjf"} ■ (42) 

R*. = (43) 

Rvx and R* are both dependent on the delay A. 

For a prescribed delay A, it is possible to find the optimum setting for the equalizer 
coefficients w op t by setting the gradient of the mean square error to zero. Applying 
5E{e 2 }/dw = 0 gives 

= R^Rvxb (44) 
This now represents the optimum w as a function of b. If this is in turn inserted into 
(40), then the following is obtained: 

E{e*} = b T • (R. - R^R^vTlW) * b - b r . C , b (45) 
This expression can be minimized. In this case, however, it is necessary to avoid the 
trivial and nonsensical solution b = 0. To ensure this, the ancillary condtoon used » 
the normalization b T -b - 1. An approach using the Lagrange multiplier then gives 

£{e 2 } = b^C-b + ACl-b^-b) (46) 

= 2C.b-+2Ab±0 (47) 
dh ^ C-b-Ab. (48) 

An eigenvalue problem thus needs to be solved. Inserting the solution into equation 
46, the following is obtained: 

E{e*} — A - b T • b — A , (49) 
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The resultant square deviation is thus equal to the eigenvalue. The solution being 
sought is thus equal to the eigenvector that belongs to the smallest eigenvalue. 

This is intended to conclude the basic illustration of this algorithm. The description of 
all variants of the algorithm would go beyond the bounds of th.s study. 

It can be assumed that all currently available ADSL modems use a variant of the 

substttute system alignment algorithm. By way of example, *e 

publication by Blade!, Moeneclaey [17] ought to form the baste « ^ the 

implementations by Alcatel, since the work for this publication was commissioned by 

Alcatel. 

3) optimization of the channel capacity: In the case of this algorithm, the channel 

capacity 

ls maximized The maximization of the channel capacity is approximately 'equivalent 
1 7e ma^izaMon of the geometric mean of the signaMo-noise ratio, (1 neglected 
Z aTu"th. togamhm,. Al. noise Influences, both externa, lane £u«d by 
o.rr brock interference, are combined and account Is a*o taken of the fact that the 
ESESSt represents recU.ngu.ar —g whfch results - -~ 
.«J. The leakaae effect becomes very obvious close to DC. There, the signe™- 
» rI«o should ideW become infinite if the NEXT function, which has a zero 
oo'inttntre T^nTsetartlng point The transfer funcUon is a,so most favoraWe 
pornt mere, u taken 8e rati0 ln the area ar0U nd DC becomes worse 

Z^Tlt ^ components are convoluted therein on account ^o the 
effect see the SNR Illustrations in Section 10). The leakage effect 'S of 
■£Z£ZX the noise, since the latter has no cyclic pref*. in contrast to the 

signal. 



alter the SNR for the individual carriers, since both the signal and the noise would be 
multiplied by the same factor, of course. 

The consideration of optimizing the channel capacity is attributed to a publication by 
Al-Dharir, Cioffi [19], where, although this optimization criterion is mentioned, only 
the geometric mean of the channel energy is maximized during implementation. In 
addition, there is again recourse to the substitute system algorithm. The central 
discovery, which now allows the channel capacity itself to be actually optimized for 
the first time, is the inclusion of the DFT s leakage effect, which was not taken into 
account previously to date. 

For the sake of simplicity, the algorithm will now be described by citing a recently 
written publication in English (the German text of the patent is somewhat 
incomprehensible on account of unagreed changes made by the patent 
department). 

-In the sequel, the steps of the iterative algorithm are described. There, New denotes 
the maximum possible number of used carriers. For baseband applications like 
DMT one has to consider conjugacy constraints to ensure that the time-domain 
signal is real. This means that the upper half of the DFT frame is dependent on the 
lower half. The DFT block length is then Noft = 2 • rW The algorithm, however, 
requires an oversampling in frequency domain by a factor of two, which leads to a 
DFT block length of Ndftc = 1 ' Ncarr- • 

The steps of the adaptation algorithm 1 : 

1 Determine the DFT transform N = (N.) of the impulse response corresponding 
' to the noise, taking into account all the spectral properties of the external 

noise. is the discrete noise power spectral density. 

2 Let Hi be proportional to the discrete channel response. The quotient HM is 
' normalized, such that E,|H,I a /£,|N,| a represents the signal-to-noise ratio SNR* 

at the input of the equalizer when using all possible carriers with the same 
transmit power 



1 <Xi) denotes the vector with components xi. 
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Ui := Hi • ^ 



3. Product Si of the channel response Hj times the equalizer response Wi In the 
DFT domain: 

Si = Hi • W, ( s2 > 

4. Determine the corresponding impulse response by applying the IFFT: 

s = IFFT(S) < 53 ) 

5. Determine the mean signal power density fraction \Nf \' of the part of the 
impulse response s that exceeds the guard interval and influences the next 
frame. 

6 Compute the product of N, times the equalizer frequency response Wji 

n£ « Ni • H * 

7 Add the noise power density components from steps 5 and 6: 

KH<+Kf (55) 

8. Compute: the autocorrelation function m that corresponds to the total noise 
power spectral density (N 4 f : 

(ni> * IFFT( <|N|| 2 > ) < 56 ) 

9. Introduce the leakage effect of the DFT by multiplying n, with a triangular 
function t\. 

X s = n { -t it i = 0,...,4/V CMT -l (57) 

According to [20], the influence of a rectangular time-domain windowing, which 
is due to the DFT operation at the receiver, corresponds to the multiplicat.on of 
the time-domain equivalent (autocorrelation function) of the noise power spectral 
density times a triangular function. In order to describe the triangular function 
approximately with the discrete time-domain vector, an oversampling in 



fluency domrin by (at M a factor of 2 is chosen. This 

Sim relies on a fluency scale wtth at least the hatf of the DMT earner 

spacing. 

10. Transform of X Into the DFT domain: 

A = FFT(X) KOO) 

11 Determine the signaUo-noise ratio SN* for every used carrier location^ No 
fancies are considered that Just resuited from ^ m ^T^L2L 
ones that are actua.ly used and are at multiples of the ong.na. earner spacing. 



SNRj = SjA it Si= |S,| 2 . i 6 U < 59 > 
U denotes the set of actually used carrier locations. 
12. Summation of all channe. capacities at all carrier locations: 

i2 = ^log 2 (l+SNR£) 

R > the palter to be ^sT^^T^Z 
multidimensional optlm.zat.on algonthm, !' ke A ^°^ c|ents . For each aeration. 
, 21) . natively modifying the ^on» c**m .en 
a„ the steps 3 to t2 have to be ^ the actual value of R is 

falls below a certain threshold, the Iteration is terminated. 

££££ the geomebic mean of all the SN*. I s U can be max^^ed. 

Batt , noise components (external noise and Inter-blocK Interference) are combined in 
Step 7 and the leakage is treated in Step 9. 

Step 5 needs to be exp, r ,d in some more 

, he detection frame «ll be posrtioned .n t,me ^ result ,ng from 

minimized. The optimum posfcon ,n ^ P ^ resulling frotn 

components of the preceding frame as „ e P distorba „ee. Four 

components of the following frame have to be consioe 



(60) 
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such exemplary impulse responses resulting from components of neighboring 
frames next to the detection frame are depicted (for illustration purposes in 
continuous ' form) in Fig. 34. The rectangular dots symbolize further impulse 
responses resulting from the following components on both sides of the detection 
frame. For every position of the impulse responses, the portion inside the detection 
frame has to be extracted by a rectangular window. The number of impulse 
responses that have to be considered is dependent on the length of the impulse 
response. Applying an FFT and afterwards computing the squared amplitudes of 
every DFT-domaln component (periodogram) yields the corresponding noise power 
spectral density (PSD). All such noise PSD contributions from the neighboring 
frames have to be summed up to obtain a mean noise PSD of the inter-block 
interference. Like for signal-independent external noise, the rectangular windowing 
of the DFT has to be considered. This means the leakage effect comes into play 
here, too. For illustration, both noise types are summed up in Step 7. before the 
leakage is incorporated in Step 9. Hereto, the detection frame has to be taken twice 
as long, like for external noise. However, the computationally more efficient solution 
would be to use a non-extended detection frame whereby the leakage is 
automatically taken into account. The less computationally efficient description has 
been preferred here, because It unveils the similarity between both noise 
components more clearly. 

Let the length of the guard interval be n G . The position of the detection frame can be 
estimated by first finding the position of a frame of length n G + 1 with the maximum 
energy content when sliding over the entire impulse response. The detection frame 
starts right after it. 

In the sequel, the procedure is again described in a more formal way. Let us first 
suppose that the time axis zero point of the total impulse response s, is located at 
the 'center of energy'. Such a time-axis shift has no influence on the results. Now, let 
the impulse response have non-zero components for i = -n v , 0, .... n n . with rv, 
being the number of precursor and n n being the number of postcursor samples. The 
mean PSD portion resulting from the part of the impulse response that 

influences the neighboring frames can be computed as follows: 

Precursors: Va u = 0 n v - a (possible positions of the impulse response, resulting 

from different samples; a. see, Fig. 34): 
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Postcursors: Vo n = 0, .... ri n - nc I a - 2: 

n n — TtQ 4- a — 2 — cr n 

As already described, we chose the frame length of the detection frame to be 4N_ 
samples long, because of similar treatment of both noise components. The real 
detection frame at a DMT receiver has. of course, only 2N«r samples. 

The mean PSD of the inter-block interference is now obtained by summing up all 
precursor and postcursor components. 



„Kn «5<'-> » FFT( S <">) , S''"' = Frr(i ( '">) • 

Upon applying the algorithm for the very small bandwidth of 32 kHz. sabiMty 
pTbtems arose during the adaptaUon. The reason for «h.s Is eerta Ny notthe 
principle of the algorithm, however, but rather probably the program ,mplemenfcrt»n 
m the case of the larger bandwidths. the parameters (sampfing rate, earner In^vaO 
were taken as a basis by ADSL. In that case, there were no P"*emsat all w* the 
equalization algorithm 3. To explain why the new equal ,zat,o^ ^^^^ 
equator settings which were not very suitable at a bandwdth of ,ust 32 kHz. a 
Zer examinaUon would have been necessary, which could not be done wrth,n the 
13 o, the study. Nevertheless, an equalizer setting was obtained w,th a sl,gh 
modification* algorilhm 3. In this case. too. an unstable response became apparent 
during adaptation, however. 



This concludes the treatment of the two central algorithms - bt allocation and time 
domain equalizer setting - in a DMT system. We had a.ready noted above that one 
problem of a DMT system is the relatively high peak-to-average ratio. The next 
section is devoted to this problem. 




neighboring 



DMT frame 
tndudlne 
guard interval 

Figure 34: Inter-block interference in neighboring frames 
3.4.3 Reduction of the peak-to-average ratio for DMT 

^Tpensation. the effecte on account of nonpar echoes are rneant to be ever 
more critical. To date, we have no experience of th, S . however. 

t mere are five basically different proposals for reducing the pealO/alues, 

, «f FFTs or even DFT transformations are of no interest h»« 
piuralrty of FFTs o _eve a|so be conceivable to have a 

application wrth ^^^'^ rang e. However, there is only one proposal 
method which operates only in tne ur i range 
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(Pat en< DeuUche TeleKom) ,n this regard, which for reasons of r^^T 
js surle ohly for 4-PSK on all earners and relatively smeli earner numbers (s 64). 

A short o«.ew of .e ^^zzs&zz rr: 

r n r^aratrrn 3s t r^ GrouP 15 * » 

somewhat more detail. 

lteratiV e peaK va.ue Clminafion tbrou* proton onto redundant subspacs 

used to eliminate the clipped signal peaks The te P proje cted. The 

onto whlct, the sfcna, profite from dipped s, g na ' ""^^ fancies. 

^ M ^T:^ZZr dpa ey » the admiss.o* values. «e««on 
l^ST* — - Nation of the ... domain -Igna. each 



time. 
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